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ABSTRACT 

We investigate the physical conditions for the growth of intermediate mass seed black holes assumed to have 
formed from remnants of the first generation of massive stars. We follow the collapse of high-cr halos with 
Tyir > 10"* K using cosmological, smooth-particle hydrodynamic (SPH) simulations in the standard ACDM 
model. During collapse of the parent halo the seed holes are incorporated through mergers into larger systems 
and accrete mass from the surrounding gas. We include a self-consistent treatment of star formation, black 
hole accretion and associated feedback processes. Even under optimistic assumptions for the seed black hole 
mass and for efficient merger rates, we find that seed holes in halos M < 10'" Mq never reach the conditions 
for critical Eddington growth. Most of the black hole growth in this regime is determined by the initial mass 
and the merger rates. Critical accretion rates are reached, albeit only after a significant delay, at the time of 
collapse (z ~ 7) for 3-4a halos of M - 10" Mq. Our results imply Mbh = 5 x 1O'^M0(M,wo/1O" Mq)"'^^ 
at the time of collapse. The required conditions of Eddington growth to explain the build-up of supermassive 
black holes 10^ Mq), as implied by Sloan quasars at z > 6, are therefore hard to meet in such a scenario. 
Without a 'jump-start' these conditions may be only achieved in extremely rare halos with Mhaio > lO'^ that 
collapsed before z ^ 6. The sub-Eddington regime in which black holes holes accrete at early time implies a 
small contribution to the reionization by miniquasar but still sufficient to cause appreciable heating of the IGM 
at z< 15- 18. 

Subject headings: quasars: general — galaxies: formation — galaxies: active — galaxies: evolution — cos- 
mology: theory — hydrodynamics 



1. INTRODUCTION 



Following the discovery o f quasars dSchmidtl 119631 : 
iGreenstein and Matthews|[l963h it was suggested that super- 
massive black holes (10^ — 10^ M©) lie at the centers of 
galaxies, and that the quasar activity is fueled by the release of 
gravitational energy from their accreted matter The remnants 
of quasar phases at early times are probably the supermassive 
black holes (SMBHs) found at the centers of galaxies in our 
local Universe. Indeed, the properties of SMBHs found at the 
centers of galaxies today are tightly coupled to those of their 
host galaxies thus providing strong observational evidence 
for a close connection between the formation and ev olution 
of ga l axies and of their central black holes ( Magorri an et alj 
[1991 iFerrarese and MerritillOOOl: iGebhardt et al.n2000l) . Re- 
markably, quasars with inferred black hole masses in excess 
of 10^ Mq have now been discovered out to z 6 (Fan et al. 
2003), posing a significant challenge for theoretical models of 
high-redshift quasar and galaxy formation. The origin, seed 
mass and the physical conditions for their subsequent growth 
to supermassive black holes remain uncertain. 

The various scenarios that have been proposed for the 
black hole 'seed' progenitors can be broadly grouped 
into two main ideas. One possible route traces the black 
hole 'seeds' to the remnants of the first generation of 
PopIII stars. Numerical simulations of the fragmentation 
of primordial clouds in the standard cold dark matter 
model suggest that these PopIII, formed of metal-free gas 
at z 20 — 30, are in d eed very massive dBromm et alj 
19991: lAbeletalJ l2000t iNakamura and Umemural I2OOII: 



to leave behind a remnant black hole seed in the range 
of 10 — IO'^Mq. Alternatively, a different mass scale of 
10^ — 10*' Mq is predicted from models in which seed 
black holes form dkectly from the collapse of low angular 
momentum, dense gas in the center of massive halos with 

virial temperatures above 10'' K ("Haehnelt and Rees| 

1993; Umemurae tail 119931: iLoeb and Rasio 199^ 



Eisenstein and Loebl 1995[: iBromm and Loebl 12003 



YoshidaetalJ 120031: iGaoetalJ l2006h and are expected 
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iKoushiappas et alJ l2004t iBegelman et all 120061) . Even 
though the disposal of angular momentum from the gas is a 
strong constraint in these models, they offer a promising sce- 
nario: seed black holes formed in this way have a 'jump-start' 
for their growth into the Sloan quasars. 

In the most commonly adopted scenario, in which the first 
black holes are traced to the first generation of stars, growing 
the seeds up to lO' Mq, (as implied by the Sloan quasars at 
z 6) requires an almost continuous accretion of gas at the 
critical, Eddington rate. If we assume that early black holes 
grow by accretion at some fraction / of the Eddington rate, 
then the black hole mass Mbh(z) = Moexp(/Af (z,Zo)/fi), 
with ts the Salpeter time and Mq and zo are the formation red- 
shift and mass, while Af (z,Zo) indicates the time between zo 
and redshift z. From this follows that with less than about 
^ 800 Myr to z = 6 and a Salpeter time of f, = 45 Myr 
(assuming a radiative efficiency of 10%), the accretion rate 
should be / > tJAt log(M(_6) /Mq) - 0.6 - 0.9 of Edding- 
ton, or equivalently an Eddington accretion rate should be 
sustained at least for / ^ 60 — 90% of the time to z = 6. 
It is an open question however whether such vigorous Ed- 
dington growth can be reached soon enough and sustained in 
the shallow potential at early times. The question of how 
much mass is accreted and therefore how much energy is 
liberated by the first miniquasars has also a significant im- 
pact for the proposed scenarios of reionizatio n and X-ray 
heating of the intergalactic medium (IGM)(e.g. iMadau et alj 
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19991 IValaseas and SiUd 1 19991 iMiralda-Escude et all 


2000; 




Venkatesan et alJ 2001 ; Wyithe and Loebl 


2003; 


Madau et alj 120041 iRicotti and Ostriked |2004< ,Zaroubi et alj 



2006h . The first black holes can potentially produce very in- 



tense ionizing radiation and due to their X-ray radiation are 
efficient at heating up the surrounding IGM. 

In this paper, we use high resolution SPH simulations of 
massive, isolated (3-4 a) halos to study the early growth and 
evolution of the first seed black holes and explore the impact 
of their evolution to reionization and X-ray heating. Semi- 
analytic models that consider plausible scenarios for the hier- 
archical ass embly and growth of massi v e black holes from 
such seeds (IVolonteri et alJl2003L l2005t iMadauet al.ll2004 
have either assumed Eddington accretion or have argued that 
quasi-spherical, 'Bondi' accretion in halos with r,,,, > 10"* K, 
where efficient cooling via hydrogen atomic lines could oc- 
cur, can sustain critical and short phases of supercritical ac- 
cretion (IVolonteri and Reell l2005h . Multi-scale simulations 
that follow the hierarchical assembly history in ACDM of the 
most massive halos forming in a 3 Gpc-* volume at z ^ 6.5 
as well as the associated black hole growth have been able 
to show such rare massive halos are goo d candidates for the 
hosts of the first quasars (iLi etalJl2006h . Here, we follow 
these earlier papers in considering the hierarchical assembly 
of SMBHs and use the approach of direct SPH cosmologi- 
cal simulations to consider the physical conditions of the gas 
inflows that feed the central black hole in its early growth. 
Thanks to the faithful tracking of the gas dynamics these 
provide a self-consistent treatment of the coupling of large- 
scale gas inflow, black hole growth and associated feedback 
processes. In particular, we will base our treatment of ac- 
cretion onto massive black hole and associa ted energy feed- 
back on the method recently developed by iDi Matteo et alJ 
(l2005h : ISpringel et alJ (l2005bl) in fuUy three-dimensional hy- 
drodynamic simulations, where, for example, we were able 
to explicitly confirm that self-regulated quasar growth can re- 
produce the observed Mbh-ct relation. We focus on under- 
standing the black hole growth history at early time and as a 
function of halo mass. We do not impose a critical Eddington 
growth phase, but instead attempt to assess the detailed black 
hole accretion history (our work can thus be seen as comple- 
mentary to Li et al.l (12006)). Our aim is to test how reliable the 
assumption of Eddington growth is likely to be, in particular 
close to the time of reionization. 

In ^ we review our simulation method and in particular 
the implementation of the black hole accretion and feedback 
model, as well as our recipe for seeding primordial minihalos 
with black holes. We then discuss the general properties and 
evolution of the black holes in simulations of halos of vari- 
ous masses (Q- We explore the implications of our results 
for BH growth scenarios as well as for the IGM reionization 
and heating in ^ In ^ we briefly summarize our main re- 
sults. Throughout, a standard ACDM cosmological model is 
assumed, with Ha = 0.7, fim = 0.3, baryon density ilh = 0.04 
and Hubble constant Hq = lOOh km/s/Mpc (h = 0.7). 

2. THE SIMULATION METHOD 

2.1. SPH cosmological simulations 

The simulations are performed with the parallel cosmo- 
logical TreePM-Smooth Particle Hydrodynamics (SPH) code 
GADGET2 ( Springe! 2005). This code has been extensively 
tested in a wide range of applications from large scale struc- 
ture formation to star formation. For the computation of 
the gravitational force, the code uses the Tree-PM method 



that combines a 'tree' algorithm for short range forces and 
a Fourier transform particle-mesh method for long range 
forces. GADGET2's implementation of SPH (e.g lMonaghanI 
fT992h uses a formulation which manifestly conserves en- 
ergy and entro py while us ing fully adaptive SPH smoothing 
lengths (iSpringel and Hernq uist 2002). Both the force com- 
putation and the time stepping of the code are fully adap- 
tiv e. Radiative coo ling and heating process are included (as 
m lKatzetalJ[T996l) . as is photoheating due to an imposed 
ionizing UV background. Detailed tests of the basic perfor- 
mance and accuracy of GADGET2 f or gravitational and hyd ro- 
dynamic problems are presented in lSpringel et alJ (l2005al) . 

The interstellar medium (ISM), star formation and super- 
novae feedback as well as black hole accretion and associated 
feedback are treated by means of sub-resolutions models. In 
particul ar, the multiphase model for sta r forming gas devel- 
oped bv lSpringel and HernquistI (l2003al) has two principal in- 
gredients: (1) a star formation prescription and (2) an effective 
equation of state (EOS). For the former we adopt a rate mo- 
tivated by observations and given by the Schmidt-Kennicutt 
Law (Kennicutt 1989), where the star formation rate is pro- 
portional to the density of cold clouds divided by the local 
dynamical time and normalized to reproduc e the star forma- 
tion r ates observed in isolated spiral galaxies dKennicutll 19891 
Il998l) . The effective EOS encapsulates the self-regulated na- 
ture of star formation due to supernovae feedback in a simple 
model for a multiphase ISM. In this model, a thermal instabil- 
ity is assumed to operate above a critical density threshold pth, 
producing a two phase medium consisting of cold clouds em- 
bedded in a tenuous gas at pressure equilibrium. Stars form 
from the cold clouds, and short-lived stars supply an energy 
of 10^' ergs to the surrounding gas as supernovae. This energy 
heats the diffuse phase of the ISM and evaporates cold clouds, 
thereby establishing a self -regulation cycle for star formation. 
Pth is determined self-consistently in the model by requiring 
that the EOS is continuous at the onset of star formation. The 
cloud evaporation process and the cooling function of the gas 
then determine the temperatures and the mass fractions of the 
two 'hot and cold' phases of the ISM, such that the EOS of the 
model can be directly computed as a function of density. The 
approach we adopt here for star formation (and the parame- 
ters we use) has akeady been shown to lead to a numerically 
converged estimate for the cosmic star formation history of 
the universe that agrees well w ith low redshift observations 
(ISpringel and Hernquisi2003bl) . 

A prescription for accretion and feedback from mas- 
sive black holes is also included (iDi Matteo et al.l 120051 : 
ISpringel et al.ll20()5bl) . Technically, we represent black holes 
by collisionless particles that grow in mass by accreting gas 
from their environments. A fraction of the radiative energy 
released by the accreted material is assumed to couple ther- 
mally to nearby gas and influence its motion and thermody- 
namic state. Our underlying assumption is that the large-scale 
feeding of galactic nuclei with gas (which is resolved in our 
simulations) is ultimately the critical process that determines 
the growth of massive black holes. Here we briefly review the 
main features and ingredients of this model. 

2.1.1. The seed black hole 

There is no clear consensus for the mass of the first black 
holes. A number of scenarios exist which suggest that rather 
than starting out as stellar mass black holes, the first black 
holes are formed as intermediate mass black holes. "Seed" 
black holes can form as the end-point of the evolution of the 
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TABLE 1 

Overview of the simulation parameters of the main runs. 
The runs are labeled with a letter indicating the type of 

RUN, and a number INDICATING THE MASS. 



label 


mass (M0//1) 




A (total, gas) 




Al 


10* 


16 


0.03 10*, 10^ 


0.05 


A2 


109 


12 


0.03 10*, 10^ 


0.05 


A3 


lO'o 


10 


0.03 10*, 10^ 


0.05 


A4 


IQii 


7.5 


0.03 5 X 10*, 5 X 10^ 


0.05 


Bl 


108 


16 


0.03 10*, 10^ 





B2 


10' 


12 


0.03 10*, 10^ 





B3 


lO'o 


10 


0.03 10*, 10^ 





B4 


10" 


7.5 


0.03 5 X 10*, 5 X 10^ 





CI 


108 


30 


0.03 10*, 10^ 


0.05 


C2 


10' 


20 


0.03 10*, 10' 


0.05 


C3 


10'° 


16 


0.03 10*, 10^ 


0.05 


C4 


10" 


9 


0.03 10*, 10^ 


0.05 



first generation of stars (e.g iB romm et al.' '1 9991: lAbel et alj 
2000; Nakamu ra and Umemur a 2001; Yos hida et al] 120031: 
Gao et al.l 120061) and in t his case a BH of mass a few 100 — 
10"^ Mp) can result (e .g iMa dau and Rees 20Qll iHe ger et al] 
120031: iRicotti and Ost riker 200^. The most massive Pop III 
stars will for m in rare high d ensity peaks of the primordial 
density field jGao et al.ll2006h . This will be the starting point 
for our study. 

We introduce collisionless 'sink' particles in the simula- 
tions to model black holes at the centers of minihalos. In 
order to achieve this we keep track of the formation of mini- 
halos by running a friends-of-friends (FOF) group finder 
For reasons of performance, and because the identity of ha- 
los changes only slowly, the group finder is not run every 
timestep but at times equally spaced in log scalefactor, with 
Alogfl = log 1.25. The group finder is run with a linking 
length of 0.16 and we select objects with a mass of 10* Mq 
and place a seed BH of Mseed = 10"^ M0 in them if they do 
not already contain a BH. In practice most of the halos of this 
mass and their black holes are formed between z = 15 and 
z = 30. 

Our choice of Mseed = 10"* M0 is on the high side for Pop III 
remnants. This mass corresponds approximately to the total 
available res ervo ir of star forming gas in these halos as found 
recently bv lGao et al.l (l2006l) . In this work, the authors could 
not accurately predict how much of this gas would end up in 
the Pop III star and although they argued that the growth prob- 
ably would stop before the stellar object reaches 10^ M0 these 
results are subject to considerable uncertainty. Hence, here 
we have explored the most optimistic case (from the point of 
view of studying early SMBH formation from stellar seeds) 
in which all of the total gas reservoir of 10"^ M0 has gone into 
the formation of first star. In ^ we will however also dis- 
cuss results for a more conservative value of Mseed = 10^ M0. 
Alternatively, these seeds can be viewed as the end results 
of direct gas collapse scenario or stellar gravitational collapse 
scenario (although these formation mechanisms are somewhat 
more speculative). 

As we shall see the black holes do not undergo much evo- 
lution at early times, so our results are rather insensitive to 
the choice of the seeding procedure. In particular different 
choices for the linking length (the over-density) give similar 



results, although choosing it too low tends to select objects 
that have not yet completed collapse and that are extended 
and non-spherical. After the seeding of the BH by this proce- 
dure, further growth of the black hole sink particles proceeds 
by gas accretion, at a rate that depends sensitively on the local 
gas conditions, or by mergers with other black hole particles. 

2.1.2. Black hole accretion & feedback model 

The accreti o n mod el is fully descri b ed an d tested in 
ISpringel et all (l2005bl) : iDi Matteo etal] (l2005h . We re- 
late the (unresolved) accretion onto the black hole to the 
large scale (resolved) gas distribution usin g a Bondi-Hoyle- 
Lyttle t on parameteriz ation (Bondi 119521: iBondi and Hovld 
Il944t iHoyle andHatleton ,1939.) . In this description, the 
accretion rate onto the black hole is given by Mb = 
47r [G^M\^p\/ {c^ -\-v\fjf/'^ where p and q are the density 
and sound speed of the ISM gas respectively, and vbh is the 
velocity of the black hole relative to the gas. However, be- 
cause of the strong dependence on sound speed of the Bondi 
accretion rate, the effective resulting from our multiphase 
ISM model will give a poor estimate of the accretion rate. If 
we make the assumption that the time spent in the cold and 
hot phase of the ISM is proportional to their respective vol- 
ume filling factors fc and //,, then formally 

Here the subscripts c and h denote quantities in the cold and 
hot phase, x is the mass fraction in the cold phase, x = fcPc/ P- 
Normally the cold phase, and hence the first term in Equa- 
tion[Tl dominates the accretion rate. We use this expression to 
calculate the accretion rates from the local conditions at the 
BH positions. We also set a limit to the accretion rate by as- 
suming the accretion rate cannot be higher than the Eddington 
rate, Mndd = (47rGMBH»ip)/(erfTc) where Wp is the proton 
mass, aj is the Thomson cross-section, and e, is the radiative 
efficiency. The latter relates the radiated luminosity, Lj, to the 
accretion rate, Mbh as — L^j (Mbhc^) , which simply gives 
the mass to energy conversion efficiency set by the amount 
of energy that can be extracted from the innermost stable or- 
bit of an accretion disk around a black hole. We will adopt 
a fixed value of = 0.1, which is the mean value found for 
a radiatively efficient (IShakura and Sunyaevll 19731) accretion 
disk on to a Schwarzschild bla ck hole. 

As in lSpringel et alJ (1200 5b'). finally, we assume that a frac- 
tion Ef of the radiated luminosity couples to the surround- 
ing gas in the form of feedback energy, viz. fifeed = Cf = 
EfEyMBHC^ ■ For simplicity, we model this energy deposition 
as thermal energy deposited isotropically in the region around 
the black hole. Lack of spatial resolution precludes us from 
considering mechanical modes of releasing the energy, like 
jets and winds. However, it is likely that any form of black 
hole feedback will lead to a fraction of the energy being ther- 
malized eventually, and that the final impact of the feedback 
depends primarily on the total amount of energy released and 
less on the form it is released in. 

2.1.3. Black Hole Mergers 
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Fig. 1. — Evolution of gas density for a representative run. Shown is tlie gas density (indicated by brigtlmess) color coded by temperature for the A3 run at 
(from top left to bottom right) redshifts of 20,15.7, 12.4 and 9.7. Width of the panels is 200 kpc (comoving). BH twice as massive as the seed are marked by red 
dots, the largest BH in each frame is enlarged and shown in black (for redshift 12 the two heaviest differ less than 20% in mass, and are both enlarged). 



When two or more objects merge to form a single dark 
matter halo, dynamical friction on their central black holes 
quickly brings them to the center of the halo where they are 
also expected to merge eventually - hierarchical black hole 
mergers contribute to the growth of the central black holes. 
However, whether black hole binaries coalesce efficiently is 
still a matter of debate. In a stellar environment, it has been 
argued that the binary hardens very s lowly (Begelman et al] 
[1980; Milo savljevic and Merrittll2003h . In gaseous environ- 
ments however binaries do coalesce mo re rapidly owing to 
strong dynamical fric tion within the gas dMakino and Funatol 
2p04; Escala et al. 2004). In our galaxy-sized simulations (e.g 
Di Matteo et al. 2005), and even more so in the cosmological 
top-hat models, it is not possible to treat in detail the problem 
of binary hardening. Because galaxies have typically large 
central concentration of gas we instead assume that two black 
hole particles merge if they come within the spatial resolution 
of the simulation (i.e. within the local SPH smoothing length) 
and their relative speed Ues below the local sound speed. This 



could introduce a dependence of the derived merger rates on 
the resolution, however we will check that this is not the case 
(see ^nA[ . 

In the final stage of the black hole mergers, the emis- 
sion of gravitational waves carries linea r momentum, which 
can cause the black ho le to recoil (iMadau et al.l 120041 : 
IVolonteri and ReesI 12006!) . Note that in our cosmological 
simulations we do not have the resolution nor the rela- 
tivistic physics required to directly calculate the ejection of 
black holes by gravitational recoil. If the recoil is larger 
than the halo escape velocity the black hole will be ejected 
from its halo (Haiman '2004t lYoo and Miralda-Escudell2004l: 
Volonteri and Rees 2006). This effect is only important for 
similar mass mergers, moreover, the proper method to cal- 
culate the recoil velocity is still uncertain, leading to several 
order of magnitude variation according to different choices of 
parameters. We will discuss the possible implications of ne- 
glecting this effect in ^ 
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2.2. Initial Conditions 

We set-up the initial conditions as isolated spherical over- 
densities that correspond to high-cr peaks in the Gaussian ran- 
dom field of cosmological density fluctuations ( Katz I99l1) . 
In this configuration, the 'top-hat' (i.e. uniform density) 
cloud is given perturbations (using the Zeldovich approxi- 
mation) with a given power spectrum P{k) oc k~^-^ appro- 
priate for the scales we study. This is the s ame tech- 
nique as that exploi t ed rec ently bv iBromm and Loeb (2003]); 
iBromm and LarsonI (|2004|) in the treatment of the primordial 
star formation problem, except that here we are interested in 
larger halos (Tyir > 10"^ K) which cool efficiently via atomic 
hydrogen lines and for which H2 cooling is not important. 

Specifically, the following procedure is employed to set 
up the particle positions and velocities: first, a regular grid 
of dark matter (DM) particles is setup. The DM particles 
are assigned perturbations using the Zeldovich approxima- 
tion, which also set s the pecuhar velocities. The amplitude 
is normalized as in Bromm et al.l (|2002) by fixing the initial 
variance cr,- of the fluctuation spectrum at the starting redshift 

Zi, 

aj^Aj^k", (2) 



such that 



1 



1 +Zvir 



1, 



(3) 



where Zvh is the target redshift of collapse (i.e. virial- 
ization). Particles from a spherical region of size R = 
(3M/47rp)'/-' are then selected and given a Hubble expan- 
sion (modified to affect a collapse at the target redshift). A 
pre-determined amount of angular momentum is imparted by 
assigning rigid rotation to the particles with a spin parame- 
ter A = L\E\^^^ /GM^/^, where L is the angular momentum, 
\E\ the total energy and M the total mass. This is necessary 
because we cannot follow self-consistently the angular mo- 
mentum evolution of the collapsing object. Gas particles are 
setup in the same way, except that they are not subject to the 
random perturbations, as on these scales pressure forces will 
have erased any initial density fluctuations. 

2.3. The simulation parameters 

Using the above procedure we can investigate different col- 
lapse scenarios by choosing the total mass M, the redshift of 
collapse Zvir and A. In Table [T] we summarize the simula- 
tion parameters of the systems we examine. They comprise 
halos of M = 10^ 10^ 10'° and lO" Mq. For Zvii- we wifl 
take the collapse redshift of a 3-(T peak (for ag — 0.9, for 
a = 0.75 these would be « 3.3-cr) of the corresponding 
mass, i.e. Zvir — 16, 12, 10 and 7.5 respectively for these ha- 
los. Note, however, that in dense regions, likely the precur- 
sors of the halos that harbor high redshift quasars, high den- 
sity peaks are over-repre sented, as found e.g. in recent nu- 
merical investigations bv lGao et al.l (l2005h (where for exam- 
ple the precursor of simulated halos reach 10^ M© at redshift 
z = 24, IO'^Mq atz = 20, 10'" Mq atz== 16 and 10" M© at 
z = 10). Thus we will also investigate a scenario with higher 
Zxdr for the above mentioned range in halo masses, correspond- 
ing to ~ 4-(T peaks. We take A — 0.03 as the standard spin 
parameter for our runs. This value is consistent with the av- 
erage spin parameter found in co smological simulations (e.g. 
iJang-Cond ell and Hernquistl2001l) . We have also run our sim- 
ulations using A = or A = 0.05 , and found no significan t dif- 
ferences in our results (see also lBromm and Loebll2003l) . For 




Fig. 2. — Projection of the gas distribution (top) and rest frame UBV 
composite rendering (bottom). Shown is a region 10 kpc wide of the A3 run 
al z= 10. The UBV panel includes stellar evolution and extinction by the 
ISM. 



the runs we take = 10^ — 10^ particles and a corresponding 
spatial resolution of a few tens to a hundred parsec. 

An important parameter for our study is the value of the 
BH feedback factor. We will consider two extreme val- 
ues: ef = corresponding to the case without feedback and 
ef 0.05, the value that has been used to reproduce the ob- 
served n ormilization of th e local Mbh-c in galaxy merger sim- 
ulations (iDi Matteo et al.i r2005). All simulations are run to a 
final redshift of about z = 6. After the collapse redshift there 
is no further infall of material or merging in these models, so 
the top-hat is not a realiable tool to investigate evolution of 
structure after this time. 

3. RESULTS 

3.1. The evolution of the central host and the seed black hole 
population 

In Figure [T] we show a short sequence of projections of the 
baryon distribution (color coded by temperature) of the A3 
run, representative of the simulations presented here. The se- 
quence starts at redshift z = 20 and ends at the redshift of col- 
lapse Zvir ^10. Each panel is 200 kpc (comoving) on the side. 
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Fig. 3. — Evolution of the total stellar mass and black hole mass for the main simulations. Different panels show the results for halos of different mass (each 
panel labeled with the corresponding mass), the solid lines show the total stellar mass, the dashed Hnes show the total BH mass. Colors indicate the type of 
simulation: blue line=A (default), green line = B (no feedback) and red line = C (higher z collapse). 
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Fig. 4. — The evolution of the black hole mass function of the A3 run 
(lO'^Mo). The plot shows the number of black holes N{Mbh > M) per halo 
of the 10'" Mq A3 run as a function of mass for four different redshifts in 
the simulation. 



As expected, the over-dense region initially expands with the 
Hubble flow at a reduced rate and subsequently turns around 
(at z ^ 16 — 17) and eventually collapses. The evolution of 
the central object, as well as the formation of substructure fol- 
lows the general scenario of structure formation in a ACDM 
cosmology. Briefly after turn around, the dark matter compo- 
nent has developed significant substructure and baryons start 
falling into the deepest potentials, black holes are formed and 



star formation eventually ensues (the only black holes shown 
in Figure [T] are the most massive and those twice the initial 
seed mass). Close to the redshift of virialization, z,.,r ^ 10 the 
halo undergoes violent relaxation, with the gas settling in the 
central regions. 

Figure |2] shows a more detailed view (10 kpc proper) of 
the gas distribution and a simulated rest-frame UBVI image 
composite just before final collapse at z ^ 10. At this point, 
what looks like a heavily obscured star bursting system fed 
by remnant bridges and filaments of HI gas, is about to merge 
with two other smaller systems. Blue colors show the in- 
tense star formation going on (the star formation rate is about 
1.5Mo/yr), indicating the general youth of the system. The 
total star mass is about 5 x 10' Mq (see Fig.O. The remain- 
ing massive black holes (that have not yet merged) and their 
associate feedback energy give rise to the faint low density 
bubbles in the gas distribution (e.g.; in the lower right object 
of the bottom panel of Fig. |2]). The total mass of black holes 
at z = 10 is 6.3 x 10* M©, with the heaviest, at the center of 
the collapsing halo, accounting for 10% of that. Further evo- 
lution of this system will form a elliptical dwarf galaxy after 
the gas reservoir is depleted (at z = 6 the star formation rate 
has markedly leveled off, see Figure O. 

Within our model the formation of seed BH mirrors the evo- 
lution of substructure in the halo. The growth of the cumula- 
tive BH mass function and the evolution in the number of BHs 
are shown in Figure |4] and Figure Irrespectively. At z = 20 
the number of BH is still low, commensurate with the lack of 
structure at masses scales of 10* M© in the parent halo. How- 
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Fig. 5. — Black hole counts, top panel: total number of BH for the 10 Mq 
runs as a function of redshift and bottom panel: BH formation rate. Blue 
drawn line shows results of the run A3 (normal Zvk), while the red dashed 
line is for the C3 run (high Zvir)- 



ever as substructure forms in the simulation, the formation 
rate of BHs increases, peaking at z = 16 (as expected this is 
close to the turnaround redshift for this halo). Most BHs have 
just formed and therefore the mass function peaks at around 
the seed mass value. At a redshift z = 12 (lower left panel of 
Figure [T]i the growth rate is offset by the merger rate and the 
number of BH reaches a maximum. The BH mass function 
has grown significantly at most mass scales. 

Further evolution of the collapsing halo decreases the num- 
ber of BHs because of the merging of sub-halos (and the cor- 
responding mergers of central BHs). This is reflected by the 
more significant growth in the high mass tail of the mass func- 
tion. Finally, at the redshift of collapse the merger rate of BH 
peaks, resulting in a rapidly declining number of BH in the 
lower mass end and a fast growth of the central massive ob- 
ject. Also shown in Figure |5] are the resulting BH counts for 
the C3 run. In this case the collapse occurs at higher red- 
shift, and the evolution of the number of BHs happens at cor- 
respondingly earlier times. It is otherwise similar to the A3 
case, with the peak formation rate in this case occurring at 
z = 24. 

As shown in Figure [3] up to about the turn-around redshift 
of the halo the total black hole mass density is higher than the 
star density, but as significant cooling of the gas ensues after 
this time, star formation becomes significant and the resulting 
star mass density quickly exceeds the black hole mass after 
z ~ 12, the growth of which has virtually saturated by z ~ 10. 
This is also reproduced in the C3 run, where the cross over 
point for the star density to exceed the black hole mass is at 
commensurately higher redshift (z ^ 22). 

We can obtain an estimate of the global black mass function 
by taking our simulations to be representative for halos of that 
mass and multiplying the corresponsing number of BHs in 
each black hole mass bin by the expected number density of 
halos of that mass and summing over the mass range of our 
simulations, 

n (Mbh >M)=Y^ mNi (Mbh > M) (4) 
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Fig. 6. — Global BH mass function at different redshift. Shown is the 
total mass function as derived from our simulations and the Press-Schechter 
formalism, for redshifts from 16 to 7.5. Shaded area gives the range of un- 
certainty from Poisson errors for z = 7.5. Similar ranges of uncertainty apply 
for the other redshifts, but are omitted for clarity. 



with number densities n, estimated using the Press-Schechter 
formalism. This is shown in Figure |6l The basic picture that 
emerges from this is that low mass BH population builds up 
early, and as heavier halos collapse the range gradually ex- 
pands to higher masses, in agreement with Figure |4] The low 
end of the mass function does not fall (or at least not as much 
as in Fig. |4]i, because the high mass halos contribute compar- 
atively less to the low mass population of BH. 

In the next section we will discuss our results on the growth 
history of the massive black hole in the center of the collaps- 
ing halo as a function of halo mass. 

3.2. BH growth and accretion history 

Figure |7] shows the evolution of the black hole mass (top 
panels) and accretion rate (bottom panels) for the most mas- 
sive black hole formed at the center of the collapsed halos. 
The major BH is traced through its most massive progeni- 
tor back to the first seed, which forms at a redshift of around 
z ~ 20 for the majority of our simulations. The accretion rate 
and the black hole mass are shown for runs for which we have 
included BH feedback and for runs with with no BH feedback 
(A1-A4 and B1-B4 in Table[B. For the 10"' Mq halo, the re- 
sults for an intermediate feedback strength (ef = 0.01) are also 
shown. 

If we compare the runs of different mass halos with and 
without feedback the impact of BH feedback becomes clear: 
even with modest levels of feedback (5% of the radiative en- 
ergy in most cases, also shown is the case of 1% for the 
10'" M0 halo) the growth of the black hole is severely lim- 
ited compared to the case with zero feedback. In the models 
with no feedback the accretion rate starts off at values of the 
order or 10^'' yr~' , much the same as the models with 
feedback. However, while in the feedback runs the average 
accretion rate stays well below the critical values (by at least 
a factor of 100), in the zero feedback case the typical accretion 
rates reach the Eddington values leading to a fast growth of the 
black hole mass. Even in this case, however, the critical ac- 
cretion phase is reached only close and subsequent to the time 
of final collapse of the halos (i.e. only at z < 10); accretion is 
sub-Eddington before that. As shown by the discrete jumps in 
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Fig. 7. — Mass and accretion rate of the central black hole for simulations of different halo mass and varying feedback strength. The mass of the parent halo 
is labeled in each panel. The top panels show the total mass as a function of redshift for the most massive black hole. Bottom panels show the corresponding 
accretion rates (drawn lines) as well as the Eddington accretion rates (dashed). In all panels, the blue line indicates simulations with a feedback strength of 
£f = 0.05 (A runs), the green dotted lines those with tf = (B runs). An intermediate case for ef = 0.01 is also shown (purple thin line) for the 10'" Mq 
simulation. The dash-dotted line in the lO" Mq panel shows results for a reduced seed mass (Mseed = 10^ Mq). 



the black hole mass evolution in Figure |7J the growth of the 
central black hole is dominated by black hole mergers up to 
the time it reaches the critical accretion phase. In models with 
BH feedback, gas accretion also becomes increasingly more 
important as the halo mass increases (and hence the central 
gas densities increase). As the accretion rate increases and the 
BH becomes heavy enough, the feedback energy associated 
with accretion, is able to heat the gas sufficiently to "shut off 
any further growth. This occurs rather abruptly. The greater 
the mass of the halo, the longer the simulation with feedback 
follows (approximately) the growth curve of the simulation 
without feedback. Hence the difference is smallest for the 
10" M0 run. Figure |7] shows that the difference in evolution 
of black hole accretion rate and hence final black hole masses 
in models with and without BH feedback increases for smaller 
halo mass runs. For the lO" Mq run the accretion rates are 
similar up until the time close to the collapse (at z ~ 7.5). 
At this point, the growth barely reaches the critical Edding- 
ton value, while lower mass runs with feedback do not reach 
Eddington growth at all. So this mass may be considered the 
minimum mass for any efficient black hole growth relevant to 
quasar activity. Still, only a fraction of the final black hole 
mass (^ 10%) results from mass accretion in this case. At 
lower mass scales the growth is more inefficient still, occur- 
ring mostly due to BH-mergers (for the 10^ Mq case less than 
1% is due to accretion). This means that the masses derived 
are sensitively dependend on the assumed seed masses. 

Note that for the runs without feedback the BH starts to 
grow at the Eddington rate eventually, not from the onset of 
BH formation, but only after a significant number of merg- 



ers have occurred. When growth at the Eddington rate does 
start, an exponential runaway growth phase occurs. In this 
case < 5% of the final mass of the BH consists of the orig- 
inal seed masses. At the end of the simulations, we see that 
the final BH masses are proportionally larger for larger halo 
masses, except for the 10" Mq simulation, because in this 
case the time after the Zvir is comparatively much shorter than 
for the other simulations, leaving too little time for all the BH 
mergers and much of the significant accretion phase to occur 
(as we have stopped all simulations at Zend = 6). 

In order to better appreciate the cause of the differences in 
the BH accretion history for the cases with and without feed- 
back, we show in Figure[8]the evolution of the physical quan- 
tities that determine the growth of the black hole in our model. 
The gas density around the BH in the model with feedback 
has lower density, as well as higher sound speed than the in 
run without feedback. This is expected due to the local heat- 
ing associated with the BH feedback which eventuall y ener- 
gizes t he gas sufficiently to drive a slow wind ( Spring el et alJ 
l2005al) . Figure [8] also shows that the relative BH-gas veloci- 
ties are higher in the feedback case: if no feedback is included 
the gas settles in a quiet equilibrium which is most efficient at 
feeding the central BH (note that the decrease in density and 
sound speed at the end of the B 1 run is a result of gas depletion 
due to runaway BH growth). The effect of BH feedback on its 
surroundings is also visible in the much stronger variations 
in density and sound speed, indicating that the ISM around 
it is being stirred, and definitely not in equilibrium. This is 
consistent with the gas density projection in Figure |2] which 
shows the faint bubbles blown in the ISM as a result of the 
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Fig. 8. — Evolution of the ISM properties around the BH. Plotted are 
(from top to bottom) BH accretion rate, local gas density, sound speed and 
the relative velocity of the BH with respect to the gas for the A3 (blue lines) 
and B3 (green dotted lines) runs. 
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Fig. 9.— Resolution test for the IO'^Mq A3 run. Plotted is the mass of 
the most massive black hole in the simulation, for runs of different resolution, 
blue dashed line: (standard) N = 10*, red drawn line: N = lO'' and green 
dotted: ^ = 5 X lO". 
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3.2.1. Resolution tests 



In order to assess possible resolution effects on our results, 
we have carried out resolution tests. Figure|9]shows the black 
hole mass for the 1O'°M0 halo with three different particle 
numbers (A^ — 10^,10^, 5 x 10^) spanning a factor of 50. The 
results show consistent black hole growth. In particular, the 
final BH mass is rather insensitive to the resolution used. The 
lowest resolution run ends up with 40% lower BH mass due 
to the fact that there is some loss in resolution of the sub- 
structure that would form extra BH at high resolution, but the 
two higher resolution runs show convergence for the final BH 
mass within 10%. Even though at high redshifts there is some 
disparity between the runs at different particle numbers, there 
is no systematic trend with particle number. Some differences 
are expected due to the stochastic nature of the BH merger 
growth in our simulations, which tend to magnify small differ- 
ences in BH trajectories at different resolution. We consider 
our results to be reliable at the medium resolution, although 
for highest mass halos we use the higher particle count. 

4. DISCUSSION 

Our results caution against the commonly made assumption 
for the early black hole growth, namely that they accrete at the 
Eddington rate (e.g. Haima n 2004 ; Yoo and Miralda-Escudl 
[2004; Madau et al. 2004; Volonteri a nd Rees 2006j). We have 
shown that this is likely an overestimate of the growth rate, 
and may consequently lead to overly optimistic mass esti- 
mates for the central black holes masses of the first quasars. 
In our simulations, seed black holes typically accrete well be- 
low the critical Eddington rates. For halos M < 10'" M0 BHs 
never reach an Eddington growth phase, or if they do, it is 
only at the time close to the collapse redshift and under the 
rather extreme assumption that any kind of black hole feed- 
back is switched off. The growth reaches the Eddington rate 
at late times only for halos of M > 10'' Mq. The final black 
hole masses at the time of collapse are proportional to the total 
halo mass with 

Mbh^Sx 10''Mo(M,,,,„/10" Mq^'' (5) 
(Fig. [Toll. For the highest mass halo this implies a black hole 
mass of Mbh ~ 10^ Mq, a couple of orders of magnitude be- 
low the values deduced for the z = 6 Sloan quasars. Our re- 
sults therefore imply that the first black holes, remnant from 
the first generation of stars, are unlikely to grow into the first 
quasars unless the Sloan quasars are hosted in halos rarer than 
3-4 a peaks. If we compare the number density in halos in- 
ferred from the space density of quasars at z = 6 (Fan et al. 
2004) this requires that their host masses be larger than those 
studied here, and of the order of 10'^ typical of an ex- 
tremely rare 5.5-(j (or higher for WMAP3) peak (note how- 
ever that this assumption corresponds to a duty cycle of order 
unity for the Sloan quasar). An extrapolation of Figure[TOl(us- 
ing Eq. |5]l gives a BH mass of ~ 10** M© for a 10'^ Mq halo 
(independent of Zvir), which falls short by a factor of > 1 with 
respect to the inferred SDSS quasar masses. If taken at face 
value. Equation |5] implies a minimum halo mass of IO'^Mq, 
which represents a > S-cr peak at z — 6. Such peaks are much 
too rare. An earlier collapse redshift for the 10'^ M0 halo 
may allow the BH mass to grow to the SDSS values (assum- 
ing the Eddington growth is sustained after the virialization of 
the parent halo). The collapse would then be required to oc- 
cur at least Az « 0.5 earlier. This shifts the 5.5-cr peak to even 
higher <t and again makes the parent halo a factor 20 — 50 too 
rare. 

These rather stringent conclusions on the inefficiency of 
seed growth are reached despite the fact that we have made 



10 



Pelupessy et al. 



rather optimistic assumptions regarding BH growth. Namely, 
we have assumed that all halos of a certain mass will host a 
seed BH (thus without taking into account possible feedback 
effects from the formation events of these BHs on nearby ha- 
los). Furthermore, the seed BHs we assume to be formed are 
massive, certainly in the upper range of masses considered 
for BH seeds from remnant Pop 111 stars. However, the rather 
slow black hole growth implied by our results makes the prob- 
lem a rather strong function of the initial choice of seed mass 
(as most of the growth occurs at the Bondi rate or via merg- 
ers instead of an exponential Eddington phase). Smaller seed 
masses than the value chosen here would imply even slower 
growth. This is illustrated in Figure|7] where , for the 10' ' Mq 
run we have also plotted a rerun of the A4 simulation, but with 
a seed mass of Mseed =10"' M0. In this case the final mass is 
more than a factor 10 smaller, as expected. Note however that 
the BH in this case also reaches Eddington growth. This again 
shows that a halo mass of 10" Mq is sufficient for Edding- 
ton growth. Larger seed masses would allow faster growth. 
However, larger seeds (e.g. 10^ M0) than the ones we have 
chosen would merely emphasize the need of a jump-start sce- 
nario for the format ion of the z = 6 Sloan quasars. Using 
detailed simulations. iLi et alJ (l2006l) have shown that halos of 
10'^ M0 could be the hosts of the first quasar. However, in 
their work it was assumed that BH would undergo growth at 
the Eddington rate, preceding the time of major mergers at 
z ~ 7 — 12. This isolated critical growth phase implies masses 
~ 10^ Mq upon entering major mergers. Furthermore they 
used fully assembled galaxies models as their initial condi- 
tions, thus not following the assembly of galaxies and their 
BH self-consistently. 

One feature of our modeling that stands out as being incon- 
ducive to BH growth is the inclusion of BH feedback. Even 
though the exact mechanism by which the AGN acts on its 
surrounding ISM is poorly understood, it is likely that the 
end-effect of AGN feedback is that a fraction of the accreted 
rest mass energy is thermalized in the surrounding medium, 
which is what we consider. The amount of feedback energy 
we use is a rather moderate 0.5% of the rest mass of the in- 
falling material. We have also shown that simulations with 5 x 
lower feedback strengths (see Fig.|7ll yield qualitatively sim- 
ilar results. Any amount of black hole feedback is therefore 
expected to affect the black hole growth rather dramatically 
regardless of the details of the feedback processes. 

The multiphase model we use for the star forming ISM 
takes into account that we have limited resolution to follow 
the densest phases of the ISM. Nevertheless, this introduces 
some dependence on the choice of model parameters. For ex- 
ample the critical density for star formation might be a func- 
tion of the chemical composition. Higher critical densities 
would give a higher mean density for the ISM, leading more 
rapid growth. It is unlikely, however, that the star formation 
model in these early objects is radically different, because star 
formation is locally seen to proceed in much the same way 
for a wide variety of environments, amongst which very low 
metallicity ones. 

For our choice of parameters almost all BH mass ends 
up in one central BH (Fig. |5]l. This is, to some extent, 
the result of our simplified criterion that leads to BH merg- 
ers. Clearly, different assumptions are likely to decrease 
the merger rate, even though, in general, the view is that 
the friction processes that bring two BHs into close vicin- 
ity of each other will operate efficiently (e.g. Begelman et alj 
Il980t IChatterjee e"tani2003h . Here we do not take into ac- 
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Fig. 10. — BH mass and halo mass. Plotted is the mass of the cen- 
tral most massive BH at the redshift of virialization of the parent halo as 
a function of halo mass. Blue triangles are the simulations of the "A" se- 
ries (our default 3-(t runs), while the red diamonds are the "C" simula- 
tions (high collapse redshift). The dotted line is the fitted relation Mbh = 
5xlO«M0(M,,„,„/lO"M0)O-^'*. 



count the effects of momentum carried away by gravitational 
waves. This so called "gravitational recoil" (GR) can expel 
BH from their parent halos instead of merging the m. This 
effect was examined in semi-analytical models by 'Haimanj 
(2004), Yoo and Miralda-Escude (2004) and more recently 
by iVolonteri and ReesI (120061) . Haiman found that GR will 
inhibit BH g rowth unless recoil veloc i ties ar e smaller than 
^ 65km/s. lYoo and Miralda-Escudd ( |2004|) . on the other 
hand, managed to grow BH to IO'Mq by z « 10 for slightly 
different assumptions regarding the GR strength and the es- 
cape velocities from halos, with about equal contributions 
from BH mergers and gas accretion processes. In either case, 
the authors did however make the assumption that the fuel- 
ing of the BH is sufficiently effective to support BH growth 
at the Eddington rate. Interestingly, Yoo and Miralda-Escud^ 
( 120041) have examined the effects of restricting the gas accre- 
tion to halos of a minimum velocity dispersion. They found 
that for the upper ranges of the threshold dispersion they con- 
sidered the BH growth was severely impeded, which is con- 
sistent with our findings. Such a threshold would be some- 
what ad-hoc in their model, whereas we find that indeed this 
restricted mode of accretion is more realistic due to the feed- 

b ack effects. 

IVolonteri and ReesI (1200 6') explore the conditions for the 
growth of seeds up to the SDSS quasars again assuming that 
this occurs at the Eddington rate. In particular, they con- 
sider the impact of dynamical effects of gravitational recoil. 
Their main conclusion is that the occurence of SMBH at early 
times is very selective. They argue that SDSS quasars could 
be hosted in halos of M ^ 10" Mq, provided that the radia- 
tive efficiency is not too high, as in this case they are likely 
not to suffer significantly from the negative contribution due 
to dynamical effects. Alternatively for growth in higher halo 
masses (10'^ Mq or so) the kick velocities need to be consid- 
ered at the low end. Our results however highlight the diffi- 
culty of sustaining Eddinton growth in the first place in halos 
of M 10" Mq, arguing for even more selective conditions 
that allow for rapid growth. 



Early growth of massive black holes 



11 





Fig. 1 1 . — Ionizing photon production. Plotted are in blue drawn lines: the production rates of ionizing photons for the stars, in green dotted: 0bh for 
the ionizing flux from all black holes and thin black drawn line: the total <!> = <j>t + 4>m- For the stellar ionized photons it is assumed that photons are liberated 
instantaneously at the formation redshift. 



Faced with the difficulty of growing seed BHs to the Sloan 
quasars, some authors propose that short phases of super- 
Eddi ngton growth can jump sta rt the growth of BHs. In partic- 
ular, IVolonteri and ReesI (l2005h proposed a scenario in which 
BH grow in thick and dense {n — 10"* cm^^, T = 8000 K) gas 
disks formed at very high redshift, at z 20 before the uni- 
verse is enriched and metal cooling can affect collapse and 
fragmentation. These authors argue that the BH should ac- 
crete gas at rates comparable to the Bondi estimate, which for 
these conditions far exceeds the Eddington rate. In our model, 
the implicit assumption for the ISM is that some enrichment 
has taken place as a result of the Poplll star formation, sig- 
nificantly enough that star formation is proceeding and feed- 
back from SN and stellar winds pressurizes the ISM, inhibit- 
ing the formati on of large and dense central disks as those in- 
ferred by Volo nteri and ReesI (l2005l) . We cannot therefore di- 
rectly probe the scenario of super-critical accretion proposed 
but ultimately the conclusion from that work also requires ex- 
tremely rare host for the Sloan quasars. 

In our simulations we have not treated the detailed physics 
of Hll regions of the first stars in the ISM suiTounding the first 
black holes at z 20 ( rather we have op timistically assumed 
efficient HI cooling). Ijohnson and Bromm (2006) have re- 
cently shown that the seed black hole growth may experience 
an additional early bottleneck due to significant delay between 
black hole formation and the onset of efficient accretion in the 
primordial Hll gas. This work strengthens our results, and 
in general the argument that it may be hard to grow the first 
quasars from the remnants PopIII black holes. 



4.1. Reionization and X-ray heating 

Several authors have recently investigated the potential 
impact of a primordial generation of BHs on the reioniza- 
tion process and the thermal evolution of the IGM. More 
specifically, high-z BHs could produce a substantial amount 
of ionizing photons due to their potentially higher lumi- 
nosity compared to stellar type sources (per u nit mass ac- 
creted) and likely higher escape fraction (e.g. 'Madau et alj 
1999; Valaaeas and Silk 1999; Miralda-Escude et al. 200^ 
iWyithe and Loebl 120031: IMadau et al.1 12004). In addition, 
their harder spectra would also induce reionization of Hell, 
otherwise unattainable with stellar type sources (with the 
exception of very massive Poplll stars). Finally, X-ray 
photons could produce a fairly homogeneous reionization 
(due to their large mean free path), albeit not co r nplete 
(e.g. IqB EQO TI IVenkate san et"aLl l200lt iMadau et al.1 l200l 
iRicotti and"O sti-iker 2004), and substantially heat the IGM. 
Most of these calculations however have also been carried out 
under the assumption that the first black holes are growing at 
the critical Eddington rate at early time. Here, we briefly in- 
vestigate the implication of the black hole accretion history 
we obtain from our simulations for ionization and X-ray heat- 
ing of the IGM. We will limit ourselves to some simple esti- 
mates, deferring a more detailed investigation to a future pa- 
per. 

The production rates of ionizing photons hv > 13.6 eV is 
plotted in Figure [TT]for the various halo masses we have sim- 
ulated. The rate of ionizing photons from stars (f>^, is calcu- 
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lated as 0* — ijj^Mi,, where is the star formation rate and 
= 3 X 10''° photons/ M0 is the yield, typical of a stellar 
population with salpeter IMF and a metallicity of about 1 % 
Solar In the case of metal-free, massive stars, -0* can be 
more than an order of magnitude higher (as discussed above 
our assumption is that in general some enrichment has taken 
place). Thus, the stellar production rate of ionizing photons 
should be regarded as a lower limit. We calculate the BH 
photon production rate assuming that all radiation comes out 
as ionizing photons with a mean energy of {hv) = 80 eV 
(such a hard spectr um is justified for smaller mass BH, e.g. 
iMadau et all (l2004l) ). This gives a production rate for BH of 
0BH = V'BH^BH, with -0BH = 1 .5 X lO*'"' photons/ . 

From Figure [TT] it is clear that the stellar contribution to 
the overall budget of ionizing photons is dominant, with the 
exception of z ~ 17 — 20 in the larger objects. Note some 
episodes of high luminosity occur at high redshift for the 
BH, when accretion rates peak momentarily before being 
quenched by BH feedback (such episodes of increased accre- 
tion are also visible in Fig.|7]i, but they are too short too affect 
the ionization much. It should be noted, however, that the es- 
cape fraction (i.e. the fraction of emitted ionizing photons that 
are able to escape in the IGM) for stellar photons is likely to 
be smaller than that for photons produced by BH, which might 
thus be relevant during certain epochs of structure evolution. 

We estimate the total production of ionizing photons «p/„ 
similar to our estimate of the total mass function, by multi- 
plying the cumulative production of ionizing radiation of each 
halo by their expected number density and summing over the 
mass range of our simulations, 

nph{t) = V'n/ / (jJidt (6) 
,. Jo 

with number densities n, estimated using the Press-Schechter 
formalism. Figure [12] shows the contribution of stars and BH 
to the cumulative number density of ionizing photons (nor- 
malized on the global hydrogen number density). The contri- 
bution from the stars under the above assumption for the ion- 
izing photon flux and with an escape fraction of 100% for both 
components is an order of magnitude larger at all redshifts. If 
the gas is distributed uniformely and does not recombine, we 
have a reionization redshift of z 12. This conclusion is not 
affected by the contribution to ionization from the BHs. 

One of the most interesting effects of the harder radiation by 
BHs is IGM heating due to X-rays. As X-rays are much more 
likely to be absorbed by Hel rather than HI (due to the higher 
cross section), direct photoionization of HI by X-rays is negli- 
gible. Nevertheless, the electron ejected by Hel photoioniza- 
tion is more likely to ionize HI than Hel because the former 
is much more abundant. The partitioning of a photoelectron 
energy between heating the gas and secondary ionization is 
a function of the gas ionization fraction. For values above 
10% most of the energy goes into heating. For this reason, 
although HI ionization by X-rays is not likely to exceed 30%, 
heating of the IGM up to temperatures of ^ few x 10^ K can 
be rapidly obtained jShuU and van S teenberg !l985h . This has 
interesting implications on the observability of the reioniza- 
tion process at high redshift. In fact, a powerful probe of the 
neutral IGM is 21 cm line emission associated with the spin- 
flip transition of the ground state of HI. In order to observe 
such a line in emission, the gas temperature needs to be higher 
that the CMB temperature. X-rays from microquasars pro- 
vide the most viable and uniformly distributed heating source 
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redshift 

Fig. 12. — Cumulative number of ionizing photons per hydrogen atom. 
Plotted as a function of redshift in the blue drawn line: total number of ion- 
ized photons per hydrogen atom produced by stars n^^/nn, in green dashed: 
"0.bh/"h and thin black drawn line: np;,/«H = (n^.* + 'J0.bh)/''h- Drawn 
and dotted lines give the 'contribution' of the halos (A1-A4) of different 
masses towards the averaging procedure. 
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Fig. 13. — Temperature due to X-ray heating. Drawn lines give the tem- 
peratures of the IGM if heated by the X-rays generated from BH accretion 
(Eq.|2) as a function of redshift for two different values of the heating frac- 
tion, /heal = 0.13 (appropiate for a low ionization fraction x = lO^**) and 
/heat = 0.6 (,r = 0.1). Dashed hne gives the CMB background temperature. 



available (e.g. ' Chen and Miralda-Escudell2004t iKuhlen et al] 
2006)- The temperature reached by the IGM subject to a soft 
X- ray background due to accretion onto BH can be estimated 
as (iMadau et al.ll2004 

T^15Yi(^\( .\ (7) 

VO.13 J VlOMo/MpcV 

with comoving accreted mass density pacc (Pacc — 
lOMo/Mpc^ at z ^ 18 for our results) and the fraction 
/heat of the X-ray energy going into heating the gas, depend- 
ing on the ionization frac tion (values for /heat range fro m 
/heat ~ 0.13 to /heat ~ 0.6. IShull and van Steenberg dl985h ). 
As we have discusses, pacc is typically smaller than pbh as 
most the black hole mass is built up by BH-mergers. 

Summing the total accretion rate over the 3-(t halos with 
Mhaio = 10 M0 as for the ionizing photon production es- 
timate, we see (Fig. [T3] l that heating by BH accretion be- 
comes important in our model at z< 15 — 18 (depending 
on /heat), so mewhat later than the redshifts z « 22 found 
by iMadau et al. (2004), whose pacc already approached a 
few hundreds Mq/Mpc^ at that redshift. The temperatures 
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we find for z < 15 — 18 is in the range 300 — 1000, which 
are so mewhat lower than those found by iKuhlen and Madaiil 
(l2005h . X -rays from SN remnants, shocks associated with 
structure formation and Ly-a photons may also aid the IGM 
heating process. 

5. CONCLUSIONS 

In this paper we have assessed the early growth and evolu- 
tion of BH seeds assumed to form with Mseed ^10"^ M0 (an 
optimistic value) as end product of the very first generation of 
Poplll stars or by direct collapse of gas at z ^ 20 — 30. We 
have used direct SPH simulations with a simple prescription 
for BH gas accretion and associated feedback to explore the 
early accretion history of seed black holes in halos of mass 
10*^ — 10" M0, representative of the typical 3-4 a halos col- 
lapsing at 20 < z < 6. We have focused on constraining the 
gas accretion history onto the seed black holes to compare 
with the widely adopted assumption that at these times seed 
BHs grow at the Eddington rate. The latter is required if such 
seed black holes are to explain the masses associated with 
z ^ 6 quasars (without postulating a phase of supercritical 
accretion as invoked by IVolonteri and ReesI (e.g. 2005)). In 
addition, Eddington accretion is generally assumed when es- 
timating the impact of mini-quasar to reionization and their 
role heating the IGM. We found that: 

• seed BHs inhabiting halos of < 10^' do not reach 
an Eddington rate accretion phase. The growth reaches 
the Eddington rate at late times only for halos of M > 
10" M0. Due the limited time spend in an Eddington 
growth phase it is difficult to explain the occurrence of 



the z ^ 6 SDSS quasars. 

• feedback from the BHs is very effective impeding BH 
growth. Unless feedback processes operate at these 
high redshifts differently such that black hole feedback 
return <C 0.005 of the rest-mass energy of the infalling 
material to the ISM, its effects leads to a significant de- 
lay for the onset of the Eddington growth phase. 

• the main growth of BHs in these objects occurs mainly 
through BH mergers. The final masses of the cen- 
tral most massive black hole scales as Mbh ^ 5 x 
10''(M/,fl/o/10" M0)°'^^ Because critical accretion 
rates are hard to reach in these halos, the final BH 
masses depend on the initial seed mass. 

• miniquasars are not likely to contribute significantly to 
reionization. 

• X-ray heating resulting from underfed miniquasars is 
however likely sufficient to cause significant IGM heat- 
ing by z ~ 15. 



The computations reported here were performed at the As- 
trophysics Theory Cluster at Carnegie Mellon University and 
at the Pittsburgh Supercomputer Center (PSC), a leading edge 
site in the NSF Shared Cyberinfrastructure. TDM acknowl- 
edges partial support for this project from NSF grant AST- 
0607819. 
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